Collapse of Axionic Domain Wall and Axion Emission 
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We examine the collapse of an axion domain wall bounded by an axionic string. It is found 
that the collapse proceeds quickly and axion domain walls disappear. However axions are emitted 
in the collapse and its energy density increases during radiation dominated era and contributes 
significantly to the present mass density of the universe. In particular the axion emitted from the 
wall can account for the dark matter in the universe for F a > 10 10 GeV. 
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I. INTRODUCTION 

The axion is the Nambu-Goldstone boson associated with the Peccei-Quinn symmetry breaking which was 

invented as the most natural solution to the strong CP problem j^] of QCD. The Peccei-Quinn symmetry breaking 
scale F a is stringently constrained by the consideration of accelerator experiments, stellar cooling and cosmology. The 
allowed range of F a /N ( axion window ) is between 10 10 GeV and 10 12 GeV, where the integer N is the color anomaly 
of Peccei-Quinn symmetry. 

Since the Peccei-Quinn symmetry is a global U(l) symmetry, global strings (axionic strings) are produced during 
spontaneous symmetry breaking. At this stage the potential for the axion field is flat, i.e., the axion is massless. 
However the axion has a mass at QCD scale through the instanton effect. The potential is written by 

V(A)=fini(l- (1) 

where is the pion decay constant and is its mass. Then the mass of the axion is given by m a ~ f^m^/Fa/N. This 
potential(0) has Z n symmetry and takes its minimum at A = 0, (f a /N)ir, 2(/ a /_/V)7r, . . . 2f a n. Z n discrete symmetry is 
spontaneously broken to produce two dimensional topological defects, i.e., axionic domain walls. The property of the 
axionic domain wall is characterized by N. Since N domain walls stretch out from each axionic string, the network 
of the string-wall systems is very complicated and survives long enough to dominate the universe for N > 1 axionic 
domain wall || . Therefore the N > 1 domain wall is not accepted cosmologically. 

One might expect that the density of the domain wall can be diluted in the inflationary universe. However, for 
the dilution mechanism to work, both the reheating temperature and the expansion rate during inflation should be 
smaller than F a Q], which is rather unnatural requirement for many inflation models. Therefore the domain wall 
problem in the axion model is serious one. 

In the case of N = 1, the domain wall is a disk bounded by the axionic string. The wall with the string boundary 
is no longer stable and it might collapse by the surface tension. In fact, since, as seen later, the surface tension is 
stronger than the tension due to the string for the domain wall with size much greater than the width of the wall, the 
dynamics of the wall is controlled by the wall tension. Therefore the string that bounds the wall cannot prevent the 
wall from collapsing. Typically the wall has size of horizon length (~ tqco) when it is formed. Therefore the time 
scale for collapse is about tQCD ~ 10~ 4 sec and the walls disappear quickly without overdosing the universe and hence 
the domain wall seems harmless in the case of N = 1 axion model. However it is expected that a number of axions 
are produced when the axionic walls collapse. If produced axions are non-relativistic or cold, their energy density 
decreases slowly (oc a -3 a: scale factor) compared with the radiation density (oc a~ 4 ). Hence the relative contribution 
of axions to the total density of the universe increases with a until the universe becomes matter dominated. 

In this paper we study the collapse of the axionic domain wall (N = 1) by the numerical integration of equations 
of motion for the axion field and estimate the energy density of axion field after the collapse. It is found that the 
energy carried by the walls is converted to the axions, which gives significant contribution to the mass density of the 
present universe and might account for the dark matter of the universe. 

The axion emission is also expected from axionic string before tQCD [g or annihilation of axionic domain walls 
[PS. However, since it is shown later that parallel two domain walls go through each other without annihilation, it 
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is unlikely that a large number of axions are produced by interaction of two domain walls. The massless axion can 
be produced by axionic strings and they acquire mass after the QCD phase transition, which might give significant 
contribution to the total density of the present universe. However there are two independent quantitative estimations 
by Davis ]8| and by Harari and Sikivie Q and, unfortunately, they are quite different. The contribution from axionic 
domain walls is expected to be at least comparable to that from axionic strings if the estimation by Harari and Sikivie 
is correct. 



II. DYNAMICS OF DOMAIN WALL 



The dynamics of the Peccei-Quinn scalar field <f> is described by the lagrangian: 

= d^*d^ + v s M) (2) 



where A is the coupling constant. As the universe cools down, Peccei-Quinn U(l) symmetry (<f> — > e l6 (j)) is sponta- 
neously broken and the scalar field has vacuum expectation value (\<p\) = F a . After global U(l) symmetry is broken, 
axionic strings are formed. From causality argument, about one axionic string is produced within the horizon at 
T ~ F a . Since the line energy density of the string is kitF 2 , the string density p st is about 4irF 2 t/t 3 ~ 4TrF^(F a /m p i) 2 
which is only ~ 1CP 13 — 10~ 17 of total density of the universe at the formation epoch. Defining <f> = \cf)\ exp(iA/ F a ), 
the lagrangian for axion field A is derived from eq.(^) as 

C a = d^Ad^A + fiml (l - cos A 



= d»Ad ll A + V w {A) (3) 

where the second term comes from QCD instanton effect which gives axion mass m a ~ f^m^/Fa- The potential 
V W (A) has a minimum at A = 0, 1ixF a and domain walls are produced between A = phase and A = 2irF a phase. 
More precisely, the axion mass has temperature dependence and increase as 

m a (T) ~ 0Am a (T = 0)(A QC d/T) 3 - 7 , (4) 

where Aqcd is the QCD scale ~ 200MeV. The domain walls are produced when the axion mass becomes greater 
than the expansion rate of the universe, i.e., m a {T\) ~ a(T\)/a. The axionic domain wall has size of about horizon 
length at the formation time ii(Ti) and has the axionic string on its boundary. The surface tension of the axionic 
domain wall is a ~ 16m a F^ fn). For wall with size i?, the ratio of surface energy to string energy is given by 

aR 2 16m a F 2 R m a R' 1 1 

Therefore the dynamics of the axionic wall with size greater than R* ~ l/(m a ) 3> t\ is determined by the wall tension 
and the string tension can be neglected. 

The cosmological evolution of the wall is determined by the surface tension and interaction of walls; the former 
makes the wall shrink and the latter cuts the wall into small pieces. In both cases the wall finally collapses by the 
surface tension. After collapse the domain walls disappear and the energy that the wall had is converted into kinetic 
energy of axion fields. If the energy of axion fields changes like non-relativistic particles as the universe expands, the 
axion energy density is 

Pa(t) = p w all(h) (jM) » ( fi ) 

where p a (t) and p wa ii are the energy density of the axion field and the axionic wall, respectively. For relativistic axion, 
the axion density decreases as ~ a~ 4 until axion becomes non-relativistic. Then eq.dq) is changed to 



where (E a ) is the average energy of emitted axion at t±. Assuming the mean distance and the mean radius of walls 
are at\ and (3t\, p wa ii is given by 
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p wa ii{ti) = , , ' = 16TO- 3 /3 2 m a F Q 2 — — — , (8) 



(ati) 3 y 45 J m pl 

where N is the relativistic degrees of freedom at t\ ~ lGeV. We expect that the numerical parameters a and /? 
are 0(1) from causality, although the precise values should be determined by the realistic simulation of cosmological 
formation of the wall-string system, which is beyond the scope of the present paper. 

Using the entropy density of the universe s(= 27r 2 7VT 3 /45) and its conservation ( sa 3 — const.), the number density 
of the axion can be written as 



s \Tim p iJ V m a 



(9) 



Then the present density of the axion is 



Pa = 1.06 x 10 6 cm- 3 a- W" 1 ^ ( -^-) ( ^1 ) . (10) 

\Tim pl J \ m a J 

Since N = 289/4, T x ~ 2GeV(F a /10GeV)" 018 and m a = 6.2 x 10- 4 eV(F a /lO w GeV)-\ 
The contribution of the axion to the density parameter f2 is given by 

where h is the Hubble constant in units of 100km/ sec /Mpc. Therefore the axion density is large enough to account 
for the dark matter in the universe unless the axion is ultra-relativistic when it is emitted. We estimate the (E a )/m a 
by numerical simulation in the next section. 

III. SIMULATION OF COLLAPSE 

In order to follow the motion of domain walls, we have solved the evolution equation of the axion field numerically. 
When a wall piece much smaller than the cosmological horizon is considered, the cosmic expansion can be ignored. 
Then the field equation under the Minkowski background is written as 



9^ d 
dt 2 



^7" V2 ^ = -7J7(^ + ^) ' ( 13 ) 



using the potentials in the equations (||) and (||). We have employed the staggered leapfrog method to solve the 
differential equation. The model parameters are chosen such that the width of the wall is equal to ten simulation 
meshes and the vacuum energy of the string is a hundred times larger than that of the wall. The variation of the 
numerical value in the latter condition does not alter our conclusion since only the wall tension governs the motion 
of walls as we have mentioned above. The boundary condition is periodic, which is useful to check the accuracy of 
our calculation since the total energy in the simulation box is conserved. The solution of a static infinite planar wall 
under \4>\ — v is used in the initial configuration. When the wall lies in yz— plane, it is expressed as 

A(x) 

— — = 7r + 2sin _1 (tanhm a x) . (14) 
Fa 

Thus the scale of the inverse axion mass represents the characteristic thickness of the axionic wall. The basic numerical 
technique is the same one in the previous paper. See the reference Q for more details. 

First we have confirmed that approaching walls that face in parallel each other pass through one another. Widrow 
showed this is true in the case of a toy sine-Gordonpotential [l||. We have reproduced the passing phenomenon by 
numerical simulations using the potential (|l3|). Fig.|l| shows the result. The initial condition is set so that the relative 
velocity is 0.05c, where c is light velocity and the separation between walls is 200 meshes. Notice that the size of one 
mesh is equal to F^ 1 in our simulation and we take m a = 0.lF a . As the time evolves, two walls become close and 
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go away without any crush. When the initial relative velocity is much larger, for example, equal to 0.5c, the pair 
annihilation of walls occurs neither. Therefore it is unlikely that annihilation of domain walls occurs and produces 
axions during the collision process. The possibility which is described in the reference should be unfavorable. 

On the other hand, in the case of the encounter of a wall with a wall edge, i.e., a string cuts the wall and the process 
of disintegration advances JlS} ]. Hence by repeated intercommutations, large walls are broken to small pieces. When 
their size becomes comparable with the thickness scale of the wall, they collapse and radiate energy as axions. In 
order to see how fast and effective the wall collapse is, we have performed the simulation of the evolution of a small 
wall piece. As a simple example we followed the time evolution of a disk wall surrounded by a circular string in the 
previous work |T^] . The result shows that the disk wall shrinks at the velocity of light and the energy of the wall is 
converted to axions. We can say that the wall collapse process is so rapid that the radiation of gravitational waves is 
hardly expected. 

In order to see the wall collapse with higher accuracy, we have performed the simulations of a strip wall. 
Figs.||(a)(b)(c) demonstrate the two-dimensional distribution of V w . The first one shows the initial configuration 
in which there is one strip wall of infinite height and 20 meshes length. As the time evolves, the edges of the strip 
approach each other, that is, the wall size decreases. Thus the wall collapse proceeds also at about the light velocity 
and the result of the three-dimensional simulation is confirmed. In the case of Fig.^|, the initial (j) is zero everywhere 
in the box. However, the initial motion may make the wall move periodically so that the death of wall might be 
avoided. We have found that this is not the case by performing the simulations with various initial motions. Even in 
the most extreme case where the strings attached to the wall edges go away each other with the light velocity, such a 
motion could do nothing more than the slight extension of the wall life time. 

Finally we have estimated the (E a )/m a in the simulation. The time evolution of the potential energy V w and 
kinetic energy of the axion field A is shown in Fig.^|(a). In these simulations, the wall edges are smoothly connected 
to the true vacuum region which is different from the disk wall case. Since the overestimation of the field gradient is 
removed by this smoothing, the quantitative analysis of energy distribution is enabled. After the wall collapse there 
remains only axionic waves which we identify as axions and the values of the potential and kinetic energy of the axion 
become constant. Since the axion waves oscillate like ~ e~ zEat , the ratio of the kinetic energy to potential energy is 
equal to (E a ) 2 /m 2 . The simulation shows (E a )/m a ~ 3, which means that the emitted axion is mildly relativistic and 
becomes non-relativistic soon by cooling due to the cosmic expansion. As mentioned before we check the accuracy of 
our simulation by the conservation of the total energy in the simulation box. In the simulation above the total energy 
is conserved within accuracy of 20%. 

In order to confirm the estimated value of {E a )/m a we have performed the simulation of higher resolution ( string 
width = 3 meshes, wall width = 30 meshes, wall length = 60 meshes) . The result (Fig.||(b)) shows (E a )/m a ~ 3 
which is quite consistent with the result of the lower resolution run. 

Our model parameters in the simulations above correspond to m a /F a — 0.1 which is much larger than the actual 
axion model ( m a /F a ~ 10~ 23 ). To check if the estimation for (E a )/m a depends on m a /F a , we have run the 
simulation of thicker wall ( string width = 1 meshes wall width = 50 meshes, wall length = 100 meshes ) in which 
ma/Fa — 0.02. Although the conservation of the total energy becomes as bad as 50% in this case, the result indicates 
that the value of (E a )/m a is affected only slightly, which means the generality of our estimation concerning the 
relativisity of emitted axions. 



IV. CONCLUSION 



In the previous section, it is found that the axionic wall collapses with time scale R/c and the mildly- relativistic 
axions remain after the collapse. The estimated (E a )/m a is about 3, which leads to the relic density of the axion 
given by 
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™ 2 *°- 01a ^{w%&) (15) 

Since a and (3 are expected to be O(l), the contribution of the axion to the present universe is large for F a > 10 10 GeV 
(compared with the baryon density fish 2 ~ 0.013 14|) and might account for the dark matter of the universe. 



However eq.(|l5|) should be compared with the density of the coherent oscillation of axion field (cold axion) [|l5|jl6| 
given by 

/ F \ 118 

^(cold)=6.5xlO- 3±0 M I ^ v ) . (16) 
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Therefore the axion from the wall has comparable to or higher density than the cold axion unless a /? is much 



less than 1. In deriving (15), we assume that the domain wall collapses rapidly after its formation. As shown in the 
simulations this assumption is quite reasonable for walls whose size is much smaller than the horizon. However, a 
large wall does not shrink until its size becomes smaller than the horizon, which leads to the increase of the density of 
axionic domain wall and as a result the density of the axion produced by the collapse of the wall. Therefore eq.(|T5|) 
may underestimate the actual axion density and the axion from axionic walls may be more important than the cold 
axion. 

Another important source of axions is an oscillating axionic string. Let us compare our result with the axion 
emission from the axionic strings. The density of the axion from the strings is given by 

fi„ (string) ~ (1 - 0.01)(F a /10 10 GeV) 118 , (17) 

where uncertainty of a factor of 100 is due to two different estimations by Davis || and Harari and Sikivie . Thus 
the importance of the axion from the wall depends on which estimation is correct. The difference comes from the 
different assumption for the energy spectrum of the emitted axion and it is hard to judge which assumption is better. 
In the case of the axionic domain wall the spectrum or average energy of emitted axions can be estimated by numerical 
simulations more easily since axions are mostly produced at the final stage of the collapsing wall whose size is much 
smaller than the cosmic scale. 

In conclusion, the relic density of the axion produced by the collapse of axionic walls is, at least, larger than the 
baryon density for F a > 10 10 GeV and accounts for a part or all of the dark matter in the universe. The axion from 
the axionic domain wall is more important than the cold axion and comparable to that from strings if Harari and 
Sikivie's estimation is correct. To make a more precise prediction of the relic density we have to know the numeric 
parameter a and j3, which is beyond the scope of the present work and will be studied in future work. 
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FIG. 1. The situation that two walls pass through each other is shown. The simulation box is one-dimensional and its size 
is 500 meshes. The vertical axis depicts the false vacuum energy by A field, V w which is normalized as V m (irF a ) = 0.02. The 
initial separation is 200 meshes and the initial relative velocity is 0.05c. Four figures correspond to the V w distribution at t — 0, 
t = 1500, t = 3000, and t = 5000 respectively. 
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FIG. 2. Evolution of a strip wall whose length is 20 meshes is shown. The size of the simulation box is 300 2 . Figures pick 
out 100 2 part of the whole plane. The value of V w is plotted and it is normalization as V w (irF a ) = 0.02. The first figure is the 
initial configuration, the second one is that at t = 10, and the last is drawn at t = 20. 

FIG. 3. Time evolution of the kinetic energy and the potential energy of the strip walls is shown. The solid lines depict the 
kinetic energy, ^A 2 and the dashed lines are the part of the potential energy by the domain wall, V w . a) The same case as 
Figure (2). b) string width = 3 meshes, wall width = 30 meshes, wall length = 60 meshes. 
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